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ABSTRACT 

Due to its proximity (9 Mpc) and the strongly bimodal color distribution of its spectroscopically 
well-sampled globular cluster (GC) system, the early-type galaxy NGC 3115 provides one of the best 
available tests of whether the color bimodality widely observed in GC systems generally reflects a 
true metallicity bimodality. Color bimodality has alternatively been attributed to a strongly nonlin- 
ear color-metallicity relation reflecting the influence of hot horizontal branch stars. Here we couple 
Subaru Suprime-Cam gi photometry with Keck/DEIMOS spectroscopy to accurately measure GC col- 
ors and a CaT index that measures the Ca II triplet. We find the NGC 3115 GC system to be 
unambiguously bimodal in both color and the CaT index. Using simple stellar population models, 
we show that the CaT index is essentially unaffected by variations in horizontal branch morphology 
over the range of metallicities relevant to GC systems (and is thus a robust indicator of metallic- 
ity) and confirm bimodality in the metallicity distribution. We assess the existing evidence for and 
against multiple metallicity subpopulations in early and late-type galaxies and conclude that metallic- 
ity bi/multimodality is common. We briefly discuss how this fundamental characteristic links directly 
to the star formation and assembly histories of galaxies. 



Subject headings: globular clusters: general 
(NGC 3115) 



galaxies: star clusters: general — galaxies: individual 
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1. INTRODUCTION 

The fact that most massive galaxies possess globular 
cluster (GC) systems that are bi modal in color has been 
recognized for over a decade (e.g. JZepf fc Ashmanlll993l : 
IQstrov et al.l lT993). The general presence in galaxy GC 
systems of two principal subpopulations is implicit in this 
finding. It is now well-established t hat the vast major- 
ity of GCs are ol d (> 10 Gyr; e.g.. IStrader etHI [20051: 
IPuzia et"aT]|2005D . which suggests that the subpopula- 
tions are separated predominantly by metallicity rather 
than age. 

Considerable effort has since been devoted to under- 
standing the origin of metallicity bimodality as it has 
significant implications for understanding the star for- 
mation and assembly hi stories of galaxies (see review by 
IBrodie fc StradeH [2006) . This follows from the notion 
that GC formation accompanies all major star formation 
events in a galaxy's history. Once formed, GCs remain 
as unchanging (compared to galaxy starlight) bright bea- 
cons, and are the accompanying witnesses to the devel- 
opment of cosmic structure. 

This general view has been questioned, however. 
Richtler (2006), showed that a hypothetical unimodal 
metallicity distribution function (MDF) can appear bi- 
modal in some optical colors w hen scatter is inc luded 
in the color-metallicity relation. lYoon et al.l ()2006|) sug- 
gested, based on their single stellar population (SSP) 
models, that a strong inflection appears in the cen- 
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tral part of the optical color-metallicity relations (be- 
tween [Z/H]= —1.5 and —0.5), driven by contribu- 
tions from horizontal branch (HB) stars. This inflec- 
tion can transform a unimodal MDF into two co lor 
peaks in their models. Recently, (jYoon et al.l 1201 laT Tbl) 
explored their expectations for ultraviolet/optical color 
combinations and concluded, again, that the underly- 
ing GC MDFs are generally unimodal. Recent investi- 
gations have found evidence for and against metallicity 
bimodality (lAlves-Brito et al.l [20111 iFoster et all 120101: 
iChies-Santos et al.l 120121 : iBlakeslee et al.l I2012T ) without 
a clear resolution of the issue. 

The most obvious objection to the idea of a unimodal 
MDF for GCs is that it conflicts with our knowledge of 
the best-studied GC system of all, that of the Milky Way. 
Here metallicities are tied to accurate abundance mea- 
surements for individual stars and the unambiguously 
bimodal metallicity distribution links dir ectly to distinct 
spati al and kinematical GC subsystems (Searlc & Zinn 
1978). Nonetheless, a generalization to other galaxies 
must rely on more limited information, such as metallic- 
ities derived from integrated spectroscopy, as we cannot 
obtain full chemical and phase space data for systems 
beyond the Local Group (although see lPota et alJl2~012l 
for the results of a color-kinematical survey of nearby 
GC systems). 

In principle, to resolve the bimodality question in a 
convincing way we would need to identify galaxies hav- 
ing the most clearly bimodal color distributions and de- 
rive metallicities from spectra for a significant fraction 
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of their GC subpopulations. Galaxies with unclear color 
bimodality tell us little about this particular issue, al- 
though they may be revealing about the complexities of 
galaxy assembly, as discussed below. If, in even one un- 
equivocal case, a bimodal fit to the MDF were to be 
strongly preferred over a unimodal one, and the propor- 
tion of metal-rich to metal-poor GCs were to be similar to 
the red/blue number ratio, this would provide convincing 
support for the principle of distinct GC metallicity sub- 
populations in galaxies. Similarly, if even one case with 
strong color bimodality turned out to be unambiguously 
unimodal in metallicity, this would provide a huge boost 
to the unimodal metallicity hypothesis. 

Here we present new spectroscopically-based measures 
of GC metallicity in a nearby galaxy, NGC 3115, that 
offers the best available test in an early-type galaxy of 
these distinct interpretations of color bimodality. NGC 
3 115, one of 11 nearb y galaxies in the spectroscopic study 
of lUsher et all (I2012T1 . is an SO galaxy hosting a populous 
GC system that is strongly bimodal in optical colors. 

The great advantage of NGC 3115 for exploring the 
questi on of metallicity b imodality is its proximity (~ 9.4 
Mpc; iTonrv et al.l I2001D . This is important because 
spectroscopic studies that generate reliable metallicity 
estimates for GCs are very challenging, even for the 
largest telescopes equipped with efficient spectrographs. 
This has led to the tendency for spectroscopically-derived 
metallicity estimates to be confined to the brighter subset 
of GCs in all but the closest galaxies. The brightest GCs 
are liable to display unimodal color and metallicity distri- 
butions because of the "blue tilt" phenomenon, or mass- 
metallicity relation, th at reflects self-enrichment for mas- 
sive metal-poor GCs.dStrader et all 120061: lHarris et al.l 
I2006t iMieske et all 12006( 1. The blue tilt tends to reduce 
the separation of the red and blue distributions with in- 
creasing brightness in the color-magnitude plane. 

For NGC 3115, which has a populo us GC system con - 
sisting of 546 ±80 GCs within 50 kpc (|Faifer et al.ll20Tl . 
we spectroscopically sample a significant fraction (~20%) 
of the ordinary GC population. Indeed, we obtain a re- 
liably measured metallicity index, independent of both 
color and SSP modeling considerations, down to the 
turnover of the GC luminosity function. 

Section 2 describes the new observations that allow us 
to characterize the metallicity distribution of NGC 3115 
GCs. Section 3 presents colors and measurements of the 
metallicity-sensitive Ca II triplet index (CaT) for these 
GCs, and our investigation of the dependence of the CaT 
on HB morphology. Section 4 is a discussion of these 
findings in the context of arguments in the literature for 
and against metallicity bimodality. Section 5 contains 
the summary and conclusions. 

2. OBSERVATIONS 

To study the metallicity distribution of the NGC 3115 
GC system we use the Subaru Suprime-Cam gi pho- 
tometry and the Keck/PEIMO S multiobject spectroscopy 
given by lArnold et al.l (|2011(1 but with the color correc- 
tions and revised zero points presented in lUsher et al.l 
(|2012H . These data were taken as part of the SLUGGS 
survejfj The spectroscopic program was designed to tar- 
get the region around the Ca II triplet at 8498, 8542 and 
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8662 A (CaT). To negate the effect of the strong sky 
lines in this spectra l region we employ a tec hnique de- 
scribed in detail in IFoster et all (|2010l 1201111 and fit a 
linear combination of stellar template spectra to the ob- 
served spectra while masking out those spectral regions 
containing significant sky lines. The fitted spectra are 
then continuum-normalized and the strength of the CaT 
index is measured from the normalized spectra. A Monte 
Carlo resampling technique is used to derive asymmet- 
ric confidence intervals for each measurement. Further 
details of the CaT observations and m easurement tech- 
niques are given m lUsher et al.l (I20I1 . 

3. GC METALLICITIES 

CaT has long been known to cor relate with metallic- 
ity for old stellar populati ons (e.g., lArmandroff fc Ziiml 
1988| ICenarro et all 12002') but issues have been raised 
about its usefulness as a one-to-one predictor of metal- 
licity because it may become less sensitive ("saturate") 
at higher metallicities ([Z/H]> —0.4) and because it 
may respond to changes in HB morphology due to the 
presence of Pasche n lines in the line index bandpasses 
(IFoster et alJl20loh . SSP modeling now shows that the 
flattening at high metallicities is mild, and that the effect 
of the HB is minimal (see below). 

Figure 1 is a plot of the highest S/N (>12 A ~ l ) CaT 
measurements for 71 GCs in NGC 3115 as a function of 
(g — i)o color, along with the corresponding histograms. 
A cut of (>8 A _1 ) increases the sample to 122 but does 
not significantly change the robustn ess of the Gaussian 
Mixt ure Modeling (GMM) statistics (|Muratov fc Gnedinl 
12010(1 . This galaxy exhibits the clearest color bimodality 
of any known GC system (> 99.9% confidence). Its DD 
value (a GMM measure of peak separation divided by 
peak widths) is 3.9. By comparison, bimodality "poster 
children" NGC 4594 and NGC 4472 have DD values of 
3.2 (B — R) and 3.0 [g — z), respectively. Over-plotted 
in the upper panel of Figure 1 is the color histogram 
for the entire NGC 3115 photometric sample, showing 
that the spectroscopically observed subset faithfully rep- 
resents the overall CC distribution. The bimodal nature 
of both the color and the CaT index distributions is un- 
ambiguous. 

Possible index sensitivity to the distribution of stars 
on the HB is investigated by measuring CaT in model 
spectra, using the same methods used for GC measure- 
ments. The models are based on those of Conroy and co l- 
laborators (|Conrov et al.l [200l IConrov fc Grrnnl [2010^ . 
but have new synthetic spectral libraries to allow for a- 
enhancement and explicit modeling of the CaT index. 
Two extreme cases are investigated in a 13 Gyr single 
stellar population, one with an extreme blue HB at all 
metallicities and the other with an all- red HB. Figure 2 
demonstrates conclusively that the CaT index is effec- 
tively insensitive to HB morphology; while the spectra 
themselves change with the addition of hot stars, the ef- 
fects are present in both the index bandpass and pseudo- 
continua and largely cancel out. Figure 2 also shows that 
a-enhancement does not change the shape of the relation 
between metallicity and the CaT index, as the enhanced 
models simply move along the relation defined by the 
solar-scaled models. CaT is thus an accurate measure of 
total metallicity, [Z/B], in GCs (where [Z/H] = [Fe/H] + 
0.934 x [a/Fe]). The relation between CaT and [Z/H] 
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Fig. 1. — CaT index measurements versus (g — i)o- The solid and 
dotted histograms are the spectroscopic (71 GCs) and the entire, 
background decontaminated, photometric (253 candidates) sam- 
ples respectively. The overall color distribution is well-represented 
by the spectroscopic sample. Overplotted on the histograms are 
the best-fit Gaussians from a standard mixture-modeling analysis. 
Bimodality is preferred over unimodality at the > 99.9% and 99.8% 
levels for the color and CaT index distributions respectively. 



[a/Fe] = 0.0 
■*■ [a/Fe] = 0.0, BHB 
-■- [a/Fe] = 0.3 




Fig. 2. — CaT index versus m e tallicity based on the 1 3 Gyr SSP 
models of IConrov eTaLl | |2009T) ; IConrov fc Gunnl 1 120101) . Dotted 
line/star symbols are a model with an extreme blue HB at all 
metallicities and no a-enhancement. Solid line/filled circles rep- 
resents a purely red HB population, again with solar [a/Fe]. The 
dashed line/square symbols have [a/Fe]=0.3. Although blue HB 
stars do influence the Paschen lines, their effect on the CaT index 
is negligible. Our CaT measurements are insensitive to [a/Fe]. 

is approximately linear and permits excellent metallicity 
discrimination over essentially the full range of observed 
GC metallicities (the models span -2.0 <{Z /H]< 0.0). 

Figure 3 is a plot of [Z/H], derived fro m CaT mea- 
surem ents using the CaT-[Z/H] relation of lUsher et al.l 
(2012), against (g — i)o color. Usher et al. compared data 
and SSP models for a large number of galaxies' GC sys- 
tems and defined a robust linear relation between CaT 
and [Z/H]. A GMM analysis finds the [Z/H] distribu- 
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Fig. 3. — GC color— metallicity relation. Metallicities are derived 
from o ur CaT measureme nt for 71 GCs, using the CaT— [Z/H] rela- 
tion of Usher et a . (2012). The dashed line is the color-metallicity 



model relation of Yoon ct al. (2011b) and the solid line is a two- 
component, best- fit linear color— metallicity relation that provides 
a better representation of the data. 

tion of the NGC 3115 GCs to be bi modal at the 99.8% 
confidence level. Overplotted are the lYoon etaH (|2011bf ) 
color-metallicity relation and a two-component linear fit 
to the NGC 3115 data. The broken linear fit provides a 
better representation of the data. A similar resu lt was 
found for M31 GCs, where iPeacock et "all (|201lD used 
ugriz optical photometry and spectroscopic metallicities 
and found color-metallicity relations that are all approx- 
im ately linear, without the sharp inflections predicted by 
the lYoon et"atl ()2011bD models. 

Figure 4 presents a comparison of the MDF derived 
from the CaT index with MDFs obtai ned by apply i ng the 
lYoon et al.l (|2011b[ ) model, and the lUsher et all (|2012f ) 
empirical color-[Z/H] relation (their equation 10), to the 
observed (g — i)o colors. While the Yoon et al. color- 
[Z/H] relation appears to provide an acceptable (al- 
though not optimal) fit to the [Z/H] data in Figure 3, 
it produces an MDF that is single-peaked with an ex- 
tended metal-poor tail. This differs dramatically from 
the bimodal MDFs de rived spectroscopically and from 
the lUsher et al.1 ()2012h relation. 

4. DISCUSSION 

Here we assess the cumulative evidence for and against 
GC metallicity bimodality and briefly discuss the rele- 
vance of GC metallicity subpopulations in unraveling the 
assembly histories of galaxies. 

4.1. Optical Spectroscopy 

Apart from the obvious example of our own Milky 
Way, where metallicity bimodality is irrefutable, there 
are several galaxies for which a reasonably large sam- 
ple of GC metallicities have been measured directly from 
their spectra. M49 (NGC 4472), t he brightest ga l axy in 
the Virgo cluster, was studied by iStrader et al.l (120071) 
who used the Keck/LRIS spectroscopy of iCohen et all 
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Measured (from CaT) 
Predicted by Yoon et al. (2011) 
Predicted by Usher et al. (2012) 




[Z/H] 

Fig. 4. — GC metallicity distribution functions. The solid line is 
the MDF derived from CaT measurements. It is bimodal at the 
99.8% confidence level. The da shed and do t ted line s are obtained 
from (g — i)o colors us ing; the [Y oon et al. (2011b) model color- 
[Z /H] relation, and the Usher et al. (2012) empirical color-[Z/H] 
relation respectively. 



(2003) to estimate metallicities from Lick indices for 47 
of its GCs. The resulting MDF was shown to be bi- 
modal, and the locations of the red and blue peaks were 
shown to be consistent with expectations based on the 
GC metallicity-host ga laxy mass relation (|Strader et al.l 
[2051 IPeng etaLl[2006l). 

More recently, Ulves Brito et al.l (|2011D used 112 high 
quality Kcck/DEIMOS spectra of confirmed GCs in 
the Sombrer o galaxy (NGC 45 9 4) to measure metallic- 
ities via the iBrodie fc Huchral (|1990fl method. Again, 
the MDF is clearly bimodal (KMM probability >90%). 
Moreover, they find evidence for a linear transformation 
between metallicity and (B — R) color. 

The main weakness of the NGC 4472 work is the small 
sample size, especially as a fraction of the total GC sys- 
tem, and that only the bright end of the luminosity func- 
tion is sampled. While the sample in NGC 4594 is larger, 
the spectra are generally of lesser signal-to-noise than 
those of the 71 selected objects in NGC 3115. More crit- 
ically, the NGC 3115 system is the definitive litmus test 
of metallicity bimodality because its color bimodality is 
much clearer than in either NGC 4472 or NGC 4594. 

NGC 5128 (Cen A) is another galax y with a well- 
studied GC system. iWoodlev et"aT1 (|2010D found 
bi/multimodality in the spectroscopically derived metal- 
licities and colors of the old GCs in their sample. How- 
ever, NGC 5128 is known to be the remnant of a recent 
merger and has a population of relativel y young GCs, 
presu med to have formed in that merger (|Beaslev et al.l 
2008). Although this galaxy is GC-rich and nearby, its 
complexity and special, post-merger status make it far 
from ideal as a bimodality benchmark. 

Two well-studied GC-populo us galaxies, M3 1 
(|Caldwell et al.l I20H and M87 (ICohen et all ImM) . 
have spectroscopic metallicity distributions that are 
not clearly bimodal, and for these the issue is not with 
sample size or signal-to-noise. Neither of these galaxies 
would fit our criterion for clear color bimodality. The 
peaks of the color distributions are poorly separated and 
the kurtoses are positive (M31: DD = 1.47, k = 0.23; 
M87: DD = 1.30, k = 0.19). Therefore, the lack of 
metallicity bimodality is to be expected. Moreover, only 



GCs significantly brighter than the mean were included 
in the Cohen e t al. s tudy. Interestingly though, in M87, 
iStrader et al.l (|2011l) found evidence from subsample 
kinematics and colors that multiple, old populations 
may be present, perhaps as a result of ongoing accretion 
of infalling satellite galaxies hosting GC systems with 
a range of masses and hence colors and metallicities. 
In general, the extended formation histories of central 
cluster galaxies like M87 may be expected to produce 
more complex GC systems, especially when studied in 
detail. M31 is also activ ely accreting satelli te galaxies 
and their associated GCs (jMackev et al.ll20~Toh . although 
it is not clear whether this recent accretion history is 
typical of or unusual for massive disk galaxies. 

An important point in the context of this work is that 
spectroscopically inferred MDFs, whether unimodal or 
multimodal, do not appear to be well-reproduced by ap- 
plication of the Yoon et al . models. For e x ample , in M31, 
using the photometry of iPeacock et al.l (I2010T) and the 
spectroscopic metallicities of iCaldwell et al.l (|20ll"l ). we 
find that the KS test gives a probability of 4 x 10 -8 that 
the metallicities predicted by the Yoon et al. models 
are drawn from the same distribution as as the spectro- 
scopic metallicities. A generic result from the Yoon et 
al. models is that they produce a single skewed metal- 
licity distribution with a broad peak, steeply declining 
at the metal- rich end and with a long tail to metal- 
poor clusters. lUsher et al. (2012) show that the observed 
[Z/H] GC distributions for 11 early- type galaxies are not 
as skewed or peaky as the Yoon et al. models would pre- 
dict. 

4.2. Near Infrared Colors 

Near-infrared colors, primarily reflecting the temper- 
atures of stars on the red giant branch, should be less 
sensitive to the distribution of stars on the HB. Thus a 
unimodal metallicity distribution should map into a uni- 
modal near-infrared color distribution, and the detection 
of bimodality in a near-infrared color distribution would 
strongly imply an underlying bimo dal metallicity distri- 
bution. iChies-Santos et all (|2012l ) argued that X-band 
photometry is key to a clean discrimination of metallic- 
ity, but their detailed modeling of photometric scatter 
showed how difficult it is to obtain the sample size and 
data quality needed to reve al the effect. The diffi culty is 
underscored in M87, where iKundu fc Zepl (|2007T l found 
clear GC color bim odality in I — H , ba s ed on NICMOS 
observations, while IChies-Santos et ail (|2012h failed to 
find the effect in g — K with a larger sample. 

Although NGC 139 9 shows significant GC bimodal- 
ity in optical colors, [Blakeslc e et al.l ([2012) found it 
to be unimodal in / — H . A worry is that their in- 
ferred MDF does not seem to be consistent with re- 
sults from the actual spectroscopic MDFs of other galax- 
ies, as these consistently have much more significant 
metal-poor subpopulations than the MDFs produced 
from photometry using nonlinear color-metallicity re- 
lations. Spectroscopically-derived metallicities are ur- 
gently needed for NGC 1399 GC system. Chies-Santos 
et al. found that the g — K color distribution of NGC 
4649 (M60) is best described by a bimodal distribution 
(KMM probability >95%), consistent with metallicity 
bimodality in this galaxy's GC s ystem. Using Spitzer 
observations, iSpitler et al.l (|2008l ) also found bimodal 
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optical/mid- IR (3.6 microns) color distributions in the 
GC systems of NGC 5128 and NGC 4594. 

4.3. Multiple GC Subpopulations and Galaxy Formation 

The well-established correlatio ns between GC 

color and h ost galaxy properties (iForbes et "all 1199 
lLarsen et al]l200lt IStrader et al.ll2006HPeng etldll^OO" 

must be linked directly to metallicity effects to establish 
GCs as chemodynamical tracers of the star formation 
and assembly processes that shaped galaxies over cosmic 
time. 

Several lines of evidence, including spatial and metal- 
licity distributions and specific frequencies, have led to 
suggestions that metal-poor blue GCs are formed in 
dark matter halos at the earliest ti mes and are associ- 
ated with metal-p oor galaxy halos (Stradcr et al.l [20051 : 
iForbes et~a l. 2012), while the metal-r ich red GCs closely 
trace the build up of galaxy bulges (|Forbes et al.l 120011 : 
IStrader et aLll2005l IPota et alJl2012T) . 

The observed bimodality might arise naturally in the 
hierarchical merging paradigm as a consequence of the 
fundamental relation between GC me tallicity and host 
galaxy mass (jMuratov fc Gn cdiu 2010), or it may be the 
result of a temporary truncation of GC formation at high 
redshift (IBrodie fc Straderi 12001) . 

Regardless of the eventual viability of these scenarios, 
the linkage of GC and galaxy properties can be expected 
to provide fundamental, unique constraints on galaxy for- 
mation in a hierarchical merging framework. 



5. SUMMARY AND CONCLUSIONS 

CaT indices and gi colors have been measured for 
71 GCs with high signal-to-noise spectra in NGC 3115, 



a nearby SO galaxy with a populous GC system that 
exhibits the most strongly bimodal color distribution 
of any known GC system. The CaT index has been 
shown to be effectively insensitive to HB morphology 
and a-enhancement over a range of metallicities from 
[Z/H]=— 2.0 to [Z/H]=0.0 using new stellar population 
models. CaT correlates closely with [Z /H] over the range 
of interest for GCs; it is an excellent proxy for metallic- 
ity in such systems. Color and CaT (and hence metal- 
licity) are shown to be bimodal for NGC 3115 GCS with 
high statistical significance (> 99.9% and 99.8%, respec- 
^y)- . . , . 

The suggestion of lYoon et al.1 (|2006l ) that the underly- 
ing metallicity distributions of GC systems are unimodal, 
and that the observed color bimodality is an artifact of 
strongly nonlinear color-metallicity transformations, is 
not confirmed in this galaxy. While this result does not 
necessarily apply to all early-type galaxies, it does lend 
considerable weight to the mounting evidence from di- 
rect MDFs that multiple GC metallicity subpopulations 
are common. The fact that a variety of optical colors 
and the easily-measured CaT index are excellent indica- 
tors of metallicity provides us with a key diagnostic for 
unravelling the star formation and assembly histories of 
galaxies. 

Data presented herein were obtained at the 
W. M. Keck Observatory, which is operated as a 
scientific partnership among Caltech, UC, and NASA. 
Based on data collected at Subaru (operated by NAOJ) 
via Gemini time exchange (GN-2008A-C-12). This work 
was supported by the NSF through grants AST-0808099, 
AST-0909237, and AST-1109878. 
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